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Abstract. Aim of this communication is to describe the first results of a work-in-progress 
regarding the chemical properties of gas and stars in the circumnuclear regions of nearby 
galaxies. Different techniques have been employed to estimate the abundances of chemical 
elements in the gaseous and stellar components of nuclear surroundings in different classes 
of galaxies according to the level of activity of the nucleus (normal or passive, star forming 
galaxies and AGNs). 



1. Introduction 

A recent comparative study of the features 
of the stellar component of the circumnuclear 
(r ~ 5pc) regions, carried out on three dif- 
ferent samples of galaxies, namely 1302 star- 
forming galaxies (SFG), 1996 active galaxies 
and 2000 normal galaxies (passive galaxies) 
byl RafanelU et al. (20 09)| has shown that SFGs 
are characterized by a spectral continuum dis- 
tribution dominated by O, B stars. On the con- 
trary, O and B stars are absent in active and 
normal galaxies which show a similar contin- 
uum in the optical spectral range, but a slightly 
different depth of the 4000A break, deeper in 
the spectra of normal galaxies (NG) than in 
the spectra of active galaxies. This difference 
in deepness of the 4000A break is a clear in- 
dication that there are more stars of spectral 
type A, namely tracers of a recent star forma- 
tion history, in the circumnuclear regions of ac- 
tive galactic nuclei (AGN) than in the same re- 
gions of not-active galaxies. This finding rein- 
forces the scenario according to which galac- 
tic activity is coiTelated to stellar formation in 



the proximity of the central engine of AGNs. 
Aim of this work is the characterization of the 
chemical composition of the gaseous compo- 
nent of the circumnuclear regions in both SFG 
and AGN. We have focused our attention on 
two different questions: 

- Is gas in which star formation is occurring 
in SFG rich in heavy elements (taking the 
Sun chemical composition as reference)? 

- Is the chemical composition of gas ion- 
ized by the central engine of AGN different 
from that of SFG? 

In our work, we have applied the CLOUDY 
photoionization code |FerIand et al. (199^ to 
the AGN and SFG samples on a grid of mod- 
els spanning a large area of the parameter 
space and obtained the coiTesponding observed 
emission lines. In particular, we have concen- 
trated our attention to Type 2 AGN, namely 
the Seyfert 2 galaxies in the redshift range 
[0.04,0.08], in order to avoid the broad emis- 
sion lines arising in the dense gaseous clouds 
close to the nucleus of the AGN and hidden by 
the torus in Seyfert 2 galaxies (S2G). Our anal- 
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ysis allows us to derive the features of the gas 
components in S2 and SF galaxies with density 
and opacity comparable to those of the gaseous 
components in observed galaxies. 



2. The samples of galaxies 

The three different samples of galaxy spectra 
studied in this work have been extracted from 
the SDSS DR7 spectroscopic dataset and 
classified on the basis of a set of classical 
spectroscopic diagnostic diagrams which 
exploit the flux ratios of spectral emission 
lines to determine the presence of AGN 
activity in the spectrum of the nuclear region 
of the galaxies. These samples have all been 
selected from a larger dataset composed of all 
galaxies observed spectroscopically and with 
redshift comprised in the range [0.04, 0.08], 
with signal-to-noise S /N > 5 in the spectral 
continuum at A - 5500A with the additional 
requirement of a signal-to-noise ratio S /N > 3 
in the strong emission lines [OIII]/15007, 
[OI]^6300, [OII]^3727, H/3, Ha, [NII]/16584, 
[SII]/l/l6717 6731, when observed. The selec- 
tion in redshift has been performed in order to 
minimize the effect on the spectra of the fixed 
aperture of spectroscopic fibers used for SDSS 
observations (see |Kewley et al. (200"6) 



for 

comparison) and allow the smallest possible 
contamination by stellar light emitted in the 
regions of the galaxies far outside from the 
nucleus. The first step of the selection process 
has been the separation of the sources without 
emission lines from the emission line galaxies. 
The 2000 galaxies with highest S/N ratios of 
this class have been chosen as members of 
the passive galaxies sample. The remaining 
emission line galaxies have been split in two 
subsamples, the first corresponding to the 
sources showing the signs of the presence 
of an AGN and the second containing SFG, 
Liners and composite galaxies by applying 
the diagnostic diagram based on the Oxygen 
lines flux ratios [Wona (2009)| using the fluxes 
measured by the SDSS and the empirical 
relation in the[OI]^6300/[OIII]^5007 vs 
[OII]A3727/[O III]^5007 diagram derived by 
|Vaona (2009)| The spectra of the galaxies se- 
lected as AGN-hosting sources in the previous 



step have been retrieved by the spectroscopic 
SDSS database and corrected for galactic 
reddening using the values of the extinction 
provided by the Nasa Extragalactic Database 
(NED) web-servic41 based on the maps and 
tools discussed in Schlege l et al. (1998)| Then, 
they have been reduced to the rest-frame 
using the value of the redshift as measured by 
the SDSS spectroscopic pipeline. The stellar 
continuum of the spectra of these presumptive 
Seyfert galaxies has been evaluated using 
STARLIGHT ( [Bruzual & Chariot (2003)D , 
and subtracted from the observed spectra in 
order to retain the emission and absorption 
features of the spectra. Afterthen, the spectral 
parameters of all emission lines with S /N > 5 
have been re-measured using an original 
method performing a multi-gaussian fit of 
the the emission lines developed by L. Vaona 
during his PhD thesis ( |Vaona (2009)| . The 
galaxies showing in their spectra Ha and Hfi 
line profiles broader than the [Olll]A5001 
profile have been rejected, and the remain- 
ing sample has been classified using the 
classical spectroscopic diagnostic diagrams 
known as Veilleux-Osterbrock diagrams 
IVeifleux & Osterbrock (l"987)| The last step 
of the selection process of the S2G galax- 
ies has been the extraction of the galaxies 
which satisfy the empirical relations in 
the [NII]/l6584/Ha vs [OIII]^5007/H6, 
[SII]^6717/Ha vs [OIII]/15007/H^ and 
[OI]/l6300/Ha vs [OIII]i5007/H6 dia grams 
as determined in ( Kewley et al. (2006) i for a 
sample of galaxies observed spectroscopically 
in the SDSS DR4: 
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The SFG galaxies have been extracted from 
the original sample using the Oxygen diagram 
and then applying again the appropriate rela- 
tions in the Veilleux-Osterbrock diagnostic di- 
agrams to the line fluxes as measured by the 
SDSS pipeline. To summarize, the final num- 
bers of members of the NG, SFG and S2G 
galaxy samples used for the following analysis 
are 2000, 1302 and 1996 respectively. 

3. Spectral continua 

The first step of the analysis has been the de- 
termination of the properties of the stellar pop- 
ulation in the close suiToundings of the galac- 
tic nucleus. The spectra of the galaxies belong- 
ing to each of the three samples described in 
the previous section have been corrected for 
galactic reddening using the extinctions pro- 
vided by the NED service and, similarly at 
what has been done for the selection of the 
S2G galaxies, reduced to the rest frame us- 
ing the spectroscopic redshift measured by the 
SDSS. Then, all spectra have been realigned 
to match the dispersion of 1 A/pixel in the 
spectral range [3600A, 9200A]. At this point, 
STARLIGHT has been used to produce the 
template spectra accounting for the purely stel- 
lar emission in the nuclear regions of the galax- 
ies. This code performs spectral synthesis on 
observed spectra, consisting in the decompo- 
sition of the spectrum in terms of a conve- 
nient superposition of a base of simple stel- 
lar populations of various ages and metallici- 
ties. The collection of stellar populations em- 
ployed to perform the fits of the input spec- 
tra in this work consists in a set of 92 stel- 
lar populations extracted from the Bruzual- 
Charlot library of evolutionary synthesis mod- 
els described [Bruzual & Chariot (2003)| with 
23 equally spaced distinct ages in the in- 
terval [10'', 1.3 • 10"'] years and 4 distinct 
metallicity values (Z = 4 ■ 10"\ 8 ■ 10"^ 2 ■ 
10"~,5 • 10"^). The template spectra, obtained 
by STARLIGHT, of all galaxies of the three 
samples have been normalized to the spectral 
flux measured at A = 5500A, since this region 
of the spectrum is devoid of emission and ab- 
sorption lines, thus being a good approxima- 
tion of the continuum component of the spec- 
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Fig. 1. Normalized average spectra for NGs 
(black curve), S2G galaxies (red curve) and 
SFGs galaxies (blue curve). 

trum. In figure [T] the comparison of the nor- 
malized average spectra for the three different 
families is shown (blue is used for SFG, red for 
S2G galaxies and black for NGs). While the 
average spectrum of SFGs galaxies is signifi- 
cantly higher than others in the shorter wave- 
lengths region, the behaviour of the spectra of 
S2G and NGs galaxies is very similar. The ra- 
tio of the average normalized spectra of S2G 
galaxies and NGs is shown in figure |2] in or- 
der to highlight the differences between their 
shapes. The average spectrum of S2G galax- 
ies appears to be systematically higher than the 
average spectrum of the NGs at wavelengths 
smaller than A = 5500A, thus indicating the 
presence of a small fraction of young stellar 
populations in the nuclear regions of galaxies 
that harbor a Seyfert 2 nucleus which are not 
found in the surrounding of the nuclei of no- 
emission Unes galaxies. 

4. Metallicity of the three samples of 
galaxies 

Three different methods have been used to es- 
timate the metallicity of the interstellar mat- 
ter and stellar populations in the central re- 
gions of galaxies of the three samples here con- 
sidered. For Seyfert 2 galaxies, a photoioniza- 
tion model for the gas in the narrow line re- 
gion of the AGN ( Vaona (2009 )) has been em- 
ployed, while the content of metals of central 
regions of SFs and normal galaxies has been 
evaluated using two empirical techniques, the 
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Ratio of 32 and NG speclra 



Fig. 2. Ratio of the normalized average spectra 
of NGs and S2G galaxies. 



P method ( |Pilyugin & Thuan (2005)| and the 
Mg2 index method ( [Bender et al. (1993)| l re- 
spectively. The last method, namely the em- 
pirical correlation between the Mg2 index and 
the metallicity, has been employed also to es- 
timate the metallicity of S2 and SFs galaxies 
and compared with the metallicity provided by 
the other methods. In both cases, the circum- 
nuclear stellar populations have been recon- 
structed using STARLIGHT. In the following 
sections these methods and the results of their 
application to the galaxy samples considered 
will be described. 

4.1. Photoionization models 

There are so many physical processes in- 
volved in the ionized gas emitting region of 
AGNs that building a model able to repro- 
duce the observed spectrum is a great chal- 
lenge. Following ( |Netzer (2008)| l, there are five 
main aspects to be taken into account for a cor- 
rect modeling of the ionized gas physical pro- 
cesses: photoionization and radiative recombi- 
nation, thermal balance, ionizing spectrum, gas 
chemical composition and clouds and/or fila- 
ments distribution. In the thermal balance we 
find the mechanical heating and the role of dust 
which is always present in the NLRs. We have 
excluded from our analysis the gas kinematics 
and winds and gas confinement mechanisms. 
A code such as CLOUDY [Ferland et al. (I998)| 
is able to solve numerically the ionization and 
thermal structure of a single cloud so it can 



be used to explore different parameters val- 
ues, to build more sophisticated models such as 
composite models with more clouds and/or dif- 
ferent geometries and time-dependent models. 
CLOUDY works with two defaults geometries, 
open and close respectively. The first one is 
based on simple photoionization field through 
a plane parallel slab with a defined ionized col- 
umn density, while in the second case the gas 
is distributed all around the source. The mod- 
els were made using CLOUDY code version 
06.02 [Ferland et al. (1998)1 We have assumed 
open geometry for all models. Once the ion- 
izing spectrum is defined (no-thermal power- 
law in our case) the input parameters are the 
density Ne/cm~^, the ionized column density 
Nc{H^) / cm~^ , the ionization parameter U and 
the metal abundances,Z/Z0. The ionization pa- 
rameter expresses the level of ionization of 
photoionized gas and is proportional to the ra- 
tio of the ionizing photon density over the gas 
density according to the relation: 



U 



Q(H) 
Attk^cNh 



(4) 



where c, the speed of light, is intro- 
duced to make U dimensionless. We built 
two kinds of photoionization models: single 
cloud models and composite models, the lat- 
ter being based on a combination of two 
different single-cloud models. In the litera- 
ture two fundamental approaches regarding 
composite models can be found: a dense 
cloud inside a medium with low density 
( [Binette etal . (1996) or combining separately 
two clouds ( |Komossa & Schulz (I997)| ). The 
evaluation of the model parameters on a grid 
of models is the only realistic way to compare 
a great number of observed spectra to the mod- 
els. Since our aim is to define the mean phys- 
ical parameters which characterize the NLRs 
using spectra in the visible range, a particular 
attention is put on the metallicity because this 
parameter is important in the investigation of 
the nature of gas and its origin. We have gener- 
ated a great number of models in order to fit all 
the observed lines. The well known and use- 
ful method of the line ratio diagrams employs 
few lines and so has limited separation power, 
so that if we reproduce all the measured lines 
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the line ratio diagrams are automatically sat- 
isfied by a large number of models. In these 
cases, parameters values are obtained by av- 
eraging the different models. Another impor- 
tant result obtained by this method is the ex- 
tension toward the ultraviolet (UV) and near 
infrared (NIR) wavelengths of the electromag- 
netic spectrum. Kx^ test has been employed to 
compare the models with the observed fluxes. 



4.1 .1 . One cloud models and the input 
parameters 

In order to explore the 5 -dimensional param- 
eter space we built, a set of single cloud 
models with the values for the parameters re- 
ported in Table [T] The ionizing spectra were 
assumed to be power-laws with v" in the 
range 10 ymi - SOA'ev, a cut-off at low en- 
ergies (v^^^) and a cut-off at high energies 
(v^^). This choice is in agreement with the 
results showed in (Komossa & Schulz (1997)|l 



and ( [Groves et al. (2004)), and supported by X 
band observations (see Mainieri et al. (2007)[ ). 
In order to assure consistency, the column den- 
sity values were checked and the simulation 
was interrupted when the number of transmit- 
ted ionizing photons was below 1%; otherwise, 
the simulation stops when the final tempera- 
ture reaches 40007*^. The basic model is sum- 
marized in Figure[3j the ionizing photons cross 
a slab with a fixed density, column density and 
dust to gas ratio (D/G) (the distance of the 
source is not required as a free parameter since 
the ionization parameter is fixed). Then the in- 
trinsic spectrum is calculated using CLOUDY: 
the calculation stops when either the fixed re- 
quired column density has been achieved or the 
temperature in the simulation has fallen below 
4000K. We assume that the intrinsic spectrum 
is reddened by the ISM, so when the models 
are compared to the observed spectra it is nec- 
essary to correct for extinction. 

4.1 .2. Set of chemical abundances 

The metallicity in HII regions and in star- 
forming galaxies may be determined by 
direct methods ( |Aller et al. (1984)} , measur- 




Fig. 3. Basic scheme of the photoionization 
models. 

Table 1. Input parameter values. Note that U , 
Ne and Nc(H+) values are reported in logarith- 
mic notation. 
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ing the line flux ratios, once the tempera- 
ture is determined, or through empirical but 
well tested metallicity sensitive calibrations 
( Pilyugin & Thuan (2005)| l, when it is not pos- 
sible to measure directly the abundances (high 
metallicity or low excitation cases). In the stel- 
lar photoionization case, the ionization struc- 
ture is known and so it is possible to obtain the 
total abundances by measuring the ionic abun- 
dances using the ionization correction factors 
method (see ( |Mathis & Rosa (1991)) and ref- 
erences therein). In the AGN, the situation is 
more difficult to handle since, when the ion- 
izing spectrum is a no-thermal power-law, the 
ionization structure is very complex: the partial 
ionized region is so large so that different ion- 
ization states are mixed and taking into account 
the ionic abundances making up the total abun- 
dances gets quite complicated. In this case, di- 
rect methods are not possible and the empiri- 
cal calibrations are obtained only by models. 
The choice of the sets of chemical abundances 
a critical step in the creation of the models be- 
cause they are critically linked to another fun- 
damental aspect of the simulations, i.e. the role 
of the dust and its composition. Once a ref- 
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erence sample of different abundances is as- 
sumed (usually, the solar set of values), a law 
describing the change in the chemical compo- 
sition and that will be used for the exploration 
of different metallicities, is required. The ob- 
servations of the ionized gas chemical abun- 
dances in HII regions and starburst galaxies 
show that all the elements, except nitrogen 
and helium, in first approximation are found 
in equal proportions. As usual, in the HII re- 
gions the metallicity is expressed in terms of 
oxygen abundance while in the visible spec- 
tra it is possible to measure three different ion- 
ization states and oxygen is the most abun- 
dant element after hydrogen and helium. Then, 
we can assume that all the elements scale di- 
rectly with the oxygen abundance, except for 
nitrogen and helium that must be treated sepa- 
rately. The nitrogen abundance, as well as he- 
lium, must be corrected for secondary produc- 
tion [WaCostasXEd^nd^ 

^ = .^[10-i'' + io"7+'"«(o/H)j_ 
H H 

^= 0.0737 + 0.0293^ (6) 

The composition and quantity of dust are 
also sources of uncertainties. The set of abun- 
dances needs to take into account both gas and 
dust composition in order to maintain the total 
assumed element abundances. As for the abun- 
dances of the elements, whose reference val- 
ues are assumed to be similar to the solar val- 
ues, the default values of abundance and com- 
position for dust are taken from the Galactic 
interstellar matter (ISM). Chosen the features 
of the dust, the abundances of the gas com- 
ponent are calculated by subtracting the dust 
abundances from the total. With the excep- 
tion of the noble gases, all abundances must 
be corrected for depletion into dust. This cor- 
rection corresponds to a variation in the abun- 
dances of elements in gaseous form. The de- 
termination of these corrections is a main issue 
of modeling since it is only possible to make 
assumptions based on the observed properties 
of the Galactic ISM. It is reasonable to sup- 
pose an increase of dust with metallicity be- 
cause at a metallicity increase corresponds an 
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Fig. 4. N/0 vs O abundance, black dashed 
line total abundances, filled black circle the 
Sun abundances, colored lines: pale blue 
D/G=1.00, cyan D/G=0.75, blue D/G=0.50, 
green D/G=0.25 

increment of the refracting elements. Another 
important assumption regards the dust to gas 
ratio value (D/G thereafter). By changing the 
relative internal dust content as a function of 
the metallicity, the depletion factors must also 
be changed, and a method to estimate the de- 
pletion factors with the variation of the dust 
content is presented in ( |Binette et al. (199"3)l l. 
In any case, some assumptions for the D/G 
and the variation of the dust with the metal- 
licity are necessary. It seems that the depletion 
into grains could be constant in many direc- 
tions of the Galaxy (see ( jVladilo (2002)| ), so 
we adopted the depletion factors reported in 
( [Groves et al. (200^) as reference values and a 
linear law for the scaling of the dust abundance 
with the metallicity. Using the conventional 
logarithmic notation, the relation between de- 
pletion factor and gas abundance is given by 
the formula: 

depliX) = log{XIH\„, - logiXIHU (7) 

In this work, it is assumed logiXj O) =cost 
for all chemical elements, with the exception 
of N and He. For these elements we follow 
( [Groves et al. (2004)l i and (?) by employing 
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a linear combination of the primary and 
secondary components of Nitrogen with the 
requirement of matching the adopted solar 
abundance patterns ( Asplund et al. (2005) i 



(Table |2]l. In the models, 4 D/G values were 
used in units of the ISM D/G = 0.25, 0.50, 
0.75, 1.00. For each D/G value, we defined 
seven sets of chemical abundances, using 
the solar oxygen abundance as reference: 
Z/Z0=O.5; 1; 1.5; 2; 2.5; 3; 3.5. We focused 
our attention on the high metallicity levels 
because nuclear gas and Seyfert galaxies 
very rarely present sub-solar metallicities 
( |Groves et al. (2006)| l. For each adopted D/G 
we changed the depletion factors set in order 
to retain the total abundances. The set of 
depletion factors is constant for each assumed 
D/G value with the exception of nitrogen 
because of the secondary production. The 
nitrogen depletion factor must change with 
the metallicity because, in agreement with our 
assumptions, if the dust increases with the 
metallicity there must be an excess of nitrogen 
in gas form due to the formula 14.1.21 The 
Nitrogen may be depleted in the carbonaceous 
grains: the involved functional groups are 
Amines (N-H, C-N), Amides (N-H), Nitriles 
(C,N triple bond stretch), Isocyanates (- 
N=C=0), Isothiocyanates (-N=C=S), Imines 
(R2C=N- R) and N itro groups (-N02 aliphatic, 
aromatic) |Kwok (2 007) Excluding the ammo- 
nia (NH3), which cannot survive in an intense 
ionizing field such as in the AGN case, all 
the mentioned compounds need their partner 
elements linearly increasing with the metallic- 
ity. In Table [3] the depletion factors used for 
D/G= 1.00 are reported, if this assumption 
is correct, a great increment of nitrogen gas 
abundance for Z/Z© > 1 .5 should be observed, 
suggesting that it is not necessary to invoke an 
excess of nitrogen in order to match the ob- 
served [A^//]/l6584 line intensity. In the Figure 
|4] there are the relations between log(N/0)gas 
and 12 + logiO I H)gas for each couple of D/G 
and Z/Z0 values. For the other elements the 
trend is constant at a given D/G value. It is 
important to stress that the depletion of the 
elements must be taken into consideration 
even if the nebula is without dust, such as 
in radiation-pressure dominated models of 



( [Groves et al. (2004)) l according to which 
the dust is blown away by the wind. In our 
models, we assumed only two types of dust, 
graphite and silicates; their reference abun- 
dances are [C/i/] = 1.22£' - 04 when D/G= 1 
in the graphite and [OIH] = 1.94£' - 04, 
[MglH] = 3.15£ - 05, [Siim = 2.82£ - 05, 
{FejH} - 2.10E - 05 for silicate respec- 
tively, while the distribution of their size is 
in the range [0.005,0.25] jim. We analyzed 
the abundances found in many Hll regions 
with direct methods in order to check the 
relation between oxygen abundance and the 
other elements. In (|Van Zee et al. (1998)| ), 



([Van Zee & Haynes (2006)K 



(Izotov & Thuan (1999)), ( jlzotov et al. (20 06) 
and ( jContini et al. (2002)) l a vast collection 
of chemical abundances were determined in 
186 HII regions spanning a range of radius 
in 13 spiral galaxies, 54 supergiant H 11 
regions in 50 low-metallicity blue compact 
galaxies, 67 HII regions in 21 dwarf irregular 
galaxies, 310 emission-line galaxies from 
the Data Release 3 of the SDSS and a sam- 
ple of 68 UV-selected intermediate-redshift 
(0 < z < 0.4) galaxies. Our set of models does 
not match very well with the trend seen in 
Figure |5] whereas there is a better fit with the 
samples used by (Groves et al. (2004)li 



jMouhcine & Contini (2002) 
iKennicutt et al. (2003)| but the log{N/0) 
trend is in any case overestimated when the 
abundances are greater than Zq. In order to 
compare our results with those obtained by 
the above mentioned authors, we adopted 
the same solar abundances, formalism and 
constraints. It is apparent that the UV-galaxies 
have a large spread in the log{N/0) values, 
so it could be possible that the HII regions 
and starburst galaxies have different sets of 
chemical abundances. This hypothesis could 
be tested by using a sample of gas abun- 
dances in nuclei of spiral galaxies to ascertain 
whether the adopted function is an upper limit 
of the real one. The final number of models 
models is 10390, 9782, 9022 and 8207 for 
D/G=0.25, 0.50, 0.75, 1.00 respectively. The 
models will be rejected unless the resulting 
spectra are classified as Seyfert 2 spectra in 
the ( [Veilleux & Osterbrock (1987)| ) diagnostic 
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Table 2. Adopted abundances and depletion 
factors (the depletion factor is defined as: 
depl = logiXIH\,, - log(X/HU). 



Fig. 5. the log(N/0) vs 12+log(0/H): filled red 
circles UV galaxies, blue crosses HII regions 
in spiral galaxies, green crosses HII regions in 
blue compact galaxies, blue *** HII regions in 
dwarf irregular galaxies and green *** emis- 
sion line galaxies, black line the adopted N/O 
function, the black circles show the adopted 
metallicity values Z/Z0, respectively 0.5, 1 
(filled circle), 1.5,2,2.5, 3,3.5. 



diagrams, following the Kewley's conditions 



( Kewley et al. (2006) I, and the synthetic fluxes 
are in the observed ranges, following Koski's 
diagrams ( |Koski (1978)| . 

4.1.3. Composite models 

The simple single cloud model described in 
the previous section is( only the first step 
of our modeling effort since in the ob- 
served Seyfert spectra two distinct c com- 
ponents are usually observed, yielding high 
and low ionization lines. A double compo- 
nent model with two clouds, each account- 
ing for either the high or the low ionization 
lines, can be created. We explored two kinds 
of composite models: Binette's models and 2 
clouds models. The first one is based on the 
( |Binette et al. (1996)l l work, with the only im- 
provements represented by the exploration of a 
larger region of the parameter space; the sec- 
ond model is based on a combination of two 
difi^erent and independent single cloud models 
( |Komossa & Schulz (1997)| , as described pre- 
viously. We considered all the possible combi- 
nations between two clouds with fixed power- 
law indexes and metallicities. The spectral 



el. 


log(X/h)a 


depl 






He 


-1.07 


0.00 


-1.07 




Li 


-10.95 


0.00 


-10.95 




Be 


-10.62 


0.00 


-10.62 




B 


-9.30 


0.00 


-9.30 




C 


-3.61 


-0.30 


-3.91 


-3.91 


N 


-4.22 


-0.30 


-4.52 


-4.52 





-3.34 


-0.24 


-3.58 


-3.71 


F 


-7.44 


0.00 


-7.44 




Ne 


-4.16 


0.00 


-4.16 




Na 


-5.83 


-0.60 


-6.43 


-5.96 


Mg 


-4.47 


-1.15 


-5.62 


-4.50 


Al 


-5.63 


-1.44 


-7.07 


-5.65 


Si 


-4.49 


-0.89 


-5.38 


-4.55 


P 


-6.64 


0.00 


-6.64 




c 
o 






J .Z.U 


J. LJ 


CI 


-6.50 


-0.30 


-6.80 


-6.80 


Ar 


-5.82 


0.00 


-5.82 




K 


-6.92 


0.00 


-6.92 




Ca 


-5.68 


-2.52 


-8.20 


-5.68 


Sc 


-8.95 


0.00 


-8.95 




Ti 


-7.10 


0.00 


-7.10 




V 


-8.00 


0.00 


-8.00 




Cr 


-6.36 


0.00 


-6.36 




Mn 


-6.61 


0.00 


-6.61 




Fe 


-4.55 


-1.37 


-5.92 


-4.57 


Co 


-7.08 


0.00 


-7.08 




Ni 


-5.77 


-1.40 


-7.17 


-5.79 


Cu 


-7.79 


0.00 


-7.79 




Zn 


-7.40 


0.00 


-7.40 





lines emitted by the first and second cloud 
are combined by a weighted mean, where the 
weight is given by the product of Hp luminosi- 
ties ratio and the solid angle ratio (a>[/L02) sub- 
tended by the clouds as seen from the source 
of the ionizing radiation (this ratio, called geo- 
metrical factor, is indicated by GF). The fol- 
lowing steps lead to the relation which has 
been applied to the combined fluxes. The fine 
luminosity is given by: 



L = 47Tr^Su/47T 



(8) 
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Table 3. Nitrogen depletion factors for the ra- 
tio D/G=1.00. Relative abundances of metals 
are in logarithmic notation. 



z/z© 


OIH 


N/H 




depl 


0.5 


-3.64 


-4.75 


-5.54 


-0.79 


1 


-3.34 


-4.22 


-4.52 


-0.3 


1.5 


-3.16 


-3.89 


-4.09 


-0.19 


2 


-3.04 


-3.66 


-3.8 


-0.14 


2.5 


-2.94 


-3.48 


-3.59 


-0.11 


3 


-2.86 


-3.32 


-3.41 


-0.09 


3.5 


-2.8 


-3.19 


-3.27 


-0.08 



A . 



where r is the inner radius of the nebula, (jjlAn 
is the covering factor, S the emission line in- 
tensity {erg crrT^). Then: 



L = An(fF = Anr^ScolAn 



(9) 



where F is the observed flux {erg s"' cm~^), d 
is the distance of the source. Defining the den- 
sity of ionizing photons {s~^ crrT^) as: 



= Q{H^)IAnr^ 



(10) 



where Q{H^) is the number of ionizing photons 
emitted from the source in a second and the 
ionization parameter as 



U = mnlriHC 
If there are two clouds: 

f/i/f/2 = rimlrlni ^ = f/ini/f/2n2(12) 
and, as a consequence: 

L2/L1 = 47Td^F2/4n(fFi = F2/F1 (13) 
Then: 

F2IF1 = UiniS 20)21 U2n2Sicoi (14) 

in terms of Hp this ratio can be written as: 

F2{H/3)/Fi{H/3) = GF ■ L^{H^) / L,{Hfi) (15) 

where L!^{Hp)IL!^{Hp) is the Hp luminosity ra- 
tio if 0J2 = C02- These quantities are calculated 
by CLOUDY for each model. If we want to 



Fig. 6. 2 clouds model, a)i and C02 are the sohd 
angles subtended by the clouds. 



find the emission line intensity relative to Hp 
then: 



{Fa)2 + {Fa)i 
F,o,{Hp) 



(16) 



where FtoAHp) is the total observed flux, and 
rearranging the terms we obtain: 



h = 



1 -I- 



F(HI1)2 



(17) 



^^2) and finally: 



h = 



(«2+^-GF.(/,)i 



1 -I- 



l'jhp) 



(18) 



GF 



If GF goes to zero only the second cloud is 
visible; on the contrary, if GF goes to infinity, 
only the first cloud is visible. Five values for 
GF have been used (0.25, 0.5, 1, 2 and 4), so 
that each couple of models provides five new 
synthetical spectra.The total number of models 
obtained with the constrains on D/G, a, Z/Zq 
fixed for each couple and after imposing the 
Kewley and Koski conditions, is about 9-10*. 

4.2. Comparison of models with 
observations 

We have compared all the measured lines com- 
ing from the observations with the synthetic 
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data provided by the models. Our goal is to 
reproduce as efficiently as possible the ob- 
served spectra with the aim of determining 
the most realistic set of input parameters. In 
( |01iva et al. (1999)l l, the authors proposed a 
new method for deriving abundances in the 
AGN narrow line region (NLR), consisting in a 
selection of "good"' models from a great sam- 
ple (they used 27000 models) which is capable 
of fairly reproducing the observed line spec- 
tra, and then slightly modifying the chemical 
abundances in order to determine the best esti- 
mates of the metallicity of the AGN. The very 
large number of synthetical spectra we have 
produced has been employed to determine the 
set of models yielding the most similar syn- 
thetic spectral features to each observed spec- 
trum with ax^ test. In this case, the;^^^ is given 
by: 



(19) 



where F^^^^^ - F^^ is the difference between the 
observed flux and CLOUDY flux, and cr, the 
error associated to the i-th line. The most prob- 
able set of input parameters can be therefore 
determined for each observed spectrum. The 
comparison has been carried out taking into ac- 
count the synthetical lines which are foreseen 
to be visible in the observed spectrum, since 
these lines would be detectable {S /N > 3 ratio) 
in the observed spectrum. Lines below the S /N 
value of the observed spectrum or not visible 
are excluded from the analysis. Two further pa- 
rameters considered in the models are the red- 
dening (Av) and the scale factor F{Hp). T\\sx^ 
test is performed using the measured fluxes not 
corrected by reddening while the flux models 
are reddened in order to achieve the observed 
F{Ha)IF(j]p) ratio. Since we have not corrected 
for the extinction using a fixed value, even if 
close to 3.1, many models show a significantly 
different value. Av is calculated applying the 
relation in ( jCardelh et al. (1989)| l: 



F{Hp) ± AF(H/3) (writing in compact form the 
fluxes FJF/}, {F, + AF„)IFp), FJiFp - AFp) 
and so on), have been adopted. Then the mean 
{{Fa/Fj}}) and the root mean square (cr) of the 
ratios are calculated, taking into account all the 
values within the range -cr < {FalFp) < cr, 
giving as a result a positive Av value. For each 
model up to 9 different synthetical spectra have 
been calculated. The scale factor allows us to 
compare directly the reddened model with the 
observed spectrum. Relative fluxes would have 
increased relative errors because of the prop- 
agation. The accepted model have to satisfy a 
predetermined significance level, suggested by 
the analysis of the reduced distribution and 
given by: 



X^/dof =x^/(nr-np) 



(21) 



where "dof stands for degree of freedom, "nr" 
is the number of measured lines and "np" the 
number of parameters used in the models. The 
minimal number of dof necessary in order to 
decide the acceptability level is given by the 
number of measured lines minus the number of 
parameters used in the models. In order to con- 
sider the model a good approximation of the 
observed spectrum, the x^/do f distribution is 
required to be similar to a gaussian distribution 
peaked on 1. If the peak is located at values 
<< 1 there is an overestimation in the fitting, 
presumably caused by the fact that very large 
error bars fit very well not significant set of pa- 
rameters of the model. On the other hand, if 
the peak is placed at values >> 1, the models 
cannot reproduce the observed spectra or the 
errors are too small. The significance level has 
been fixed, after several tests, between 5% and 
10% and the assumed flux errors at 2cr level. 
The best fit models are associated to the min- 
imum value, indicated with {xl,i„), with the 
assumption that if a model is accepted, all the 
models falling within the assumed confidence 
region are considered as acceptable as well. 
This condition can be expressed as: 



Av = 7.22- 



, ( F{Ha) \ ( F(Ha) \ 



(20) 



X 



A mm 



<A 



(22) 



Nine different values of F(Ha)/F{Hj3), ob- 
tained combining F{Ha) + AF(Ha) with 



where A is the value selected for a given 
probability |Molla & Hardy (2002)| (25% level 
of confidence is adopted in this work). For each 
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fitted spectrum there is a group of accepted 
models; the best values of the parameters are 
evaluated by averaging the single values over 
all the accepted models. In this analysis, 7 free 
parameters for the single cloud models (SC) 
have been considered. As a consequence, a 
minimum number of measured lines equal to 
7 and 11 to fit such models is required, while 
at least 12 measured lines and 7 parameters 
for Binette's models. About 50% of the spec- 
tra shows less than 12 measured lines, mean- 
ing that the 2C models can be compared to 
only the 50% of the original sample of galax- 
ies but with the great advantage of considering 
only the higher S/N spectra. A spectrum with 
a low number of measured lines has usually 
low S/N, i.e. then big errors on the measured 
spectral lines, so, even though it is easier to fit 
the models with different parameters, these fits 
have little or no scientific validity. 

4.2.1. The test 

In order to compare the different kinds of mod- 
els, the spectra are splitted into two groups, 
namely spectra with nr > \2 and spectra with 
nr < 12. Since the D/G ratio is not considered a 
parameter of the model, the models with differ- 
ent D/G values will be treated separately. Such 
classification leads to four families of models, 
one for each D/G value. The number of input 
parameters for each kind of model are summa- 
rized in Table m 

4.2.2. Seyfert 2 

On a total of 2153, 1 141 spectra have nr < 12 
and the remaining nr > 12. For this reason, 
the comparison between SC and 2C models 
will be done only over 1012 spectra. In Table 
|5] the number of fitted spectra for each D/G 
value is reported; in more details, the first col- 
umn contains the D/G value, second to fourth 
columns contain the fitted spectra with SC 
models and the last column the fitted spectra 
with 2C models. The total percentage of spec- 
tra which passed the test are ~ 65% and ~ 50% 
for the SC and 2C models respectively. The 
percentage of spectra successfully fitted with 









-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3 









-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Yig.l. Idof distribution for SC models, up, 
and 2C models, down, and the different D/G 
values, red 0.25, blue 0.50, green 0.75 and grey 
1. 00. 

nr < 12 is 80%, while ~ 45% is the fraction of 
spectra fitted with nr > 12. Both two 2C mod- 
els accounts for only 50% of the spectra fitted, 
since this work was carried out with the goal of 
achieving the best fitting and not the maximum 
number of spectra fitted. 

The^^/t/o/ distributions are shown in fig- 
ure |2] The distributions of the SC models are 
slightly skewed toward lower values for the 
D/G=l and D/G=0.75 models, while the 2C 
models show tails for larger values even if their 
overall shapes are more regular The significant 
difference between the models is evident com- 
paring the distribution of absolute since this 
statistics represents the goodness of the overall 
fitting procedure (the input model parameters, 
the choice of the significance level and errors) 
and it is possible to discriminate between two 
models which are able to fit the same spectrum 
by looking at the value. 

The cumulative distributions (Figure [8]l of 
X^ Idof show clearly the great difference in the 
obtained fits between the two kinds of models 
used in this work. Within each type of mod- 
els, the D/G value yielding the minimum x^ 
values is 1 even if the difference is not signi- 
ficative, while the larger overall number of fit- 
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Table 4. The input parameters of the simulations for single cloud, double clouds and Binette's 
models. 



Parameter 


SC 


2C 


Binette 


power-law index, a 


yes 


yes 


yes 


metallicity, Z/Z© 


yes 


yes 


yes 


density, Ne 


yes 


yes 


fixed 


column density, Nc(H+) 


yes 


yes 


yes 


ionization parameter, U 


yes 


yes 


yes 


density 2nd cloud, Ne2 


no 


yes 


fixed 


column density 2nd cloud, Nc(H+)2 


no 


yes 


calculated 


ionization parameter 2nd cloud, Ut 


no 


yes 


calculated 


geometrical factor, Gf 


no 


yes 


yes 


Reddening, Av 


yes 


yes 


yes 


Scale factor, F(H^) 


yes 


yes 


yes 


number of parameters 


7 


11 


7 



Table 5. Number of fitted spectra. 



D/G 


SC 


SC nr < 12 


SC nr > 12 


2C 


0.25 


1249/2153 


853/1141 


396/1012 


404/1012 


0.50 


1375/2153 


912/1141 


463/1012 


486/1012 


0.75 


1445/2153 


930/1141 


515/1012 


554/1012 


1.00 


1426/2153 


939/1141 


487/1012 


561/1012 



ted spectra is obtained with D/G= 0.75 1. 
More interesting is the comparison of be- 
tween the matched fitted spectra with single 
and 2C models. For each D/G value (starting 
from 0.25) there are respectively 289, 356, 410 
and 396 shared spectra, whose distributions 
are dramatically different. The differences can 
be appreciated looking at the cumulative distri- 
butions (Figure|9l): spectra having good S/N are 
best fitted using 2C models. Many measured 
lines and much greater number of spectra could 
be fitted increasing the parameter resolution, in 
particular for the power-law index and the ion- 
ization parameter, at the cost of a much heavier 
computational load. 

Table |6] contains the percentage of failures for 
each measured line of the observed spectra. 
The data show again that the 2C models fit 
more efficiently all the lines. Anyway, these 
results should be cautiously handled since the 
percentage of successful fits are calculated av- 
eraging over few spectra, so that a single mea- 
sure could modify drastically the result (in par- 
ticular for the [ArlV] and [FeVII] lines). The 




Fig. 8. cumulative distribution for SC mod- 
els, thin lines, and 2C models, thick lines. 



lines mostly affected by this kind of issue are 
[A^/]5200, [Fey//]572 1,6087, //e/5876 and 
[Ar///]7135. 

4.2.3. Metallicities of Seyfert 2 galaxies 

In Table |7] are summarized all the mean values 
of the parameters with their root mean square, 
cr, respectively for the single cloud and 2C 
models. The power-law index a decreases to 
-1.8 with the increasing of the D/G ratio, prob- 
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Table 6. Percentage of line failures for single cloud and double cloud models. 





SC 








2C 








Une 


0.25 


0.50 


0.75 


1.00 


0.25 


0.50 


0.75 


1.00 


[0//]3727 


6 


6 


7 


5 


3 


2 


2 


1 


[A'e///]3869 


2 


2 


1 


1 


2 








3 


rA'p///13968 


1 


1 


2 


3 


5 


2 


1 


1 


[5//]4074 


8 


18 


18 


25 



















1 


1 


1 














[0///]4363 


10 





7 


8 


3 


6 


6 


8 


HeII4686 


2 


1 


1 


3 





1 








[Ar/V]4711 


























[Ar/V]4740 


33 


25 


25 


20 














//84861 


























[0///]4959 


























[O///]5007 


2 


3 


3 


3 














[A'/] 5200 


67 


62 


59 


57 


34 


29 


32 


34 


[Fel///]5721 


50 


60 


67 


67 


14 


10 


30 


33 


HeI5876 


19 


17 


16 


15 


8 


10 


9 


11 


[Fey//]6087 


62 


61 


66 


63 


9 


15 


15 


65 


[O/]6300 


7 


6 


7 


6 


2 


3 


5 


5 


[0/]6363 


2 


1 


1 


1 








1 





[A'//]6548 


























H„6563 


























[A'//]6584 


2 


3 


4 


5 














[5//] 67 16 


2 


3 


3 


4 








1 


1 


[5//]6731 


4 


4 


4 


4 


8 


6 


4 


1 


[Ar///]7135 


7 


11 


16 


22 


3 


8 


26 


27 


[0//]7325 
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scatter; 1.5 -;- 2.5 Z/Zq is the range of the fit- 
ted values with the SC models while for the 
2C model the estimate obtained is 1.0 2.0 
Z/Zq. The high metallicity of the gas indicates 
that the AGN fuel is evolved material, which, 
in turn, supports the idea that the gas has a lo- 
cal origin in the vast majority of the observed 
Seyfert 2 galaxies. The densities are in agree- 
ment with the measured values: SC and SC 
models produce the same values of densities, 
both very similar to the measured sulfur den- 
sity. 

The column densities are very similar in 
the SC models and in the first cloud of the 
2C models, while the second cloud of the lat- 
ter model has a lower column density and a 
lower ionization parameter too, which allows 
us to treat such cloud as a ionization bounded 
case. The ionization parameter in the SC mod- 
els is almost constant for different D/G values: 




Fig. 9.x^ cumulative distribution for shared fit- 
ted spectra with the usual meaning of the sym- 
bols. 



ably because of the decreasing gas metallicity 
content (a large fraction of metals is depleted 
into grains). In fact, a low metallicity of the gas 
implies a much slower cooling and a weaker 
ionizing spectrum is necessary to obtain the 
same degree of ionization. The total metallic- 
ity is well defined in all models with negligible 
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Table 7. Parameters output 



One cloud models 


0.25 




0.50 




0.75 




1.00 




Parameter 


mean 


cr 


mean 


cr 


mean 


cr 


mean 


cr 


a 


-1.5 


0.3 


-1.6 


0.3 


-1.7 


0.3 


-1.8 


0.3 




2.0 


0.6 


1.9 


0.6 


1.9 


0.5 


1.8 


0.5 


loe(Ne) 


2.3 


0.7 


2.3 


0.7 


2.2 


0.6 


2.1 


0.7 


log{Nc{H+)) 


20.5 


0.7 


20.3 


0.6 


20.0 


0.5 


19.9 


0.4 


U 


-3.1 


0.2 


-3.1 


0.2 


-3.0 


0.2 


-3.0 


0.3 


Av 


1.4 


0.6 


1.5 


0.6 


1.5 


0.6 


1.6 


0.6 


2C models 




a 


-1.3 


0.3 


-1.4 


0.3 


-1.6 


0.3 


-1.8 


0.3 




1.6 


0.5 


1.7 


0.4 


1.7 


0.4 


1.7 


0.4 


logiNe) 


2.2 


0.8 


2.2 


0.8 


2.1 


0.7 


2.2 


0.8 


!og(Nc(H+)) 


20.2 


0.4 


20.2 


0.4 


20.2 


0.4 


20.2 


0.4 


U 


-2.7 


0.4 


-2.6 


0.4 


-2.6 


0.4 


-2.5 


0.4 


log{Ne)2 


3.0 


0.9 


2.8 


0.9 


2.7 


0.9 


2.5 


0.8 


log(Nc(H+))2 


19.7 


0.4 


19.7 


0.3 


19.7 


0.3 


19.7 


0.3 


U2 


-3.3 


0.2 


-3.3 


0.2 


-3.3 


0.2 


-3.3 


0.3 


Av 


1.2 


0.4 


1.2 


0.4 


1.3 


0.4 


1.4 


0.4 



log(U) - -3. In the 2C the high ionization 
cloud shows a little variation in the ionization 
parameter with respect to D/G values, the mean 
value being logiU) - -2.6, while a very simi- 
lar value is found for the low ionization cloud 
(log(U) = -3.3). In general, the amount of ion- 
ization in the SC models is intermediate be- 
tween the values found in the first and second 
clouds of the 2C models. Another result, which 
can be explained with an increasing dust con- 
tent, is the increase of the absorption Av with 
D/G. There is an average difference of 0.2 mag 
between the SC and 2C models. In Figure [TOl 
the values of the Av found by models for all 
the fitted spectra and obtained with all the mod- 
els are shown against the Av values estimated 
by Balmer decrement. The Av models appear 
to be slightly larger then the estimated values 
for both SC and 2C models. This effect can 
be explained by the fact that the Ha I Hp ratio 
obtained by CLOUDY, in the vast majority of 
cases, is significantly lower than the canonical 
value of 3.1. 

The temperatures obtained by the models are 
shown in Table [8] The [OH] and [SH] tem- 
peratures values are consistent with those mea- 
sured by R02 and RS2t ratios, while the [OHI] 
temperature is lower than the measured one 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 5.5 4.0 
Av mjig (Halmer deci^ment Ha|Hb=3.1) 

Fig. 10. Comparison between the observed val- 
ues of the absorption Av and the values pro- 
duced by CLOUDY. 

but nonetheless consistent with the photoion- 
ization case. 

4.3. Metallicity for star-forming galaxies 

The determination of chemical abundances in 
the HII regions of SFs galaxies is compli- 
cated by the difficulty of reliable measures 
of the temperature of the HII clouds. The 
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Table 8. Model temperatures for the observed and simulated spectral lines in single cloud and 
double clouds models. 





0.25 




0.50 




0.75 




1.00 




SC 


Te 


cr 


Te 


(T 


Te 


tr 


Te 


tr 


[OIII] 


8274 


1421 


8522 


1437 


8879 


1387 


9374 


1379 


[0/7] 


7813 


1145 


8142 


1184 


8568 


1130 


9136 


1192 


[Sin 


7026 


1038 


7282 


1163 


7608 


1220 


8021 


1527 


2C 




[OIII] 


10923 


3257 


11172 


3189 


11225 


2937 


11535 


2853 


[oin 


10020 


3353 


10390 


3169 


10361 


2806 


10736 


2658 


[Sin 


9219 


3715 


9560 


3567 


9334 


3294 


9759 


3170 


[Olllh 


9282 


1118 


9298 


1085 


9516 


1103 


9679 


1140 


[Oil], 


8788 


1078 


8906 


1016 


9222 


1020 


9464 


1006 


[SII]2 


7995 


1418 


8049 


1353 


8307 


1450 


8450 


1480 



auroral lines, such as [OIII]4363, [SIII]6312 
and [NII]5755, used to evaluate the temper- 
ature are usually weak and in low-excitation 
metal-rich HII regions they are often too weak 
to be detected. For this reason, the direct 
method is almost useless in the vast majority 
of the cases. The paper by (Pagel et al. (1979) i 
suggested that some empirical calibrations, 
based on collection of observations, between 
the nebular line fluxes could be used to de- 
rive directly the oxygen abundance or the 
temperature of the HII regions. In more 
details, the parameter defined as R23 - 

([O//]3727.[O///]4959,5007) ^^^^ ^^^.^^^j^ 

measure the Oxygen abundances. This ratio 
was recalibrated using photoionization mod- 
els and confirmed through other observations. 
Nonetheless, the comparison between the re- 
sults obtained from different samples of galax- 
ies remains very difficult because of the high 
heterogeneity of the model parameters and 
data. In all cases, the one-dimensional calibra- 
tions are systematically wrong as showed in 
( [Pilyugin & Thuan (2005)l l. The only R23 in- 
dex does not remove all the degenerations, so a 
more general two-dimensional parametric cal- 
ibration, called the P-method, was proposed 
( Pilyugin & Thuan (2005)) l, introducing a new 
parameter, called the excitation parameter and 
defined as: 

P,3- ([O///]4959,5007) 

([O//]3727h-[O///]4595,5007) ^ ^ 



From a collection of high precision 
measurements of 700 HII regions, in 
(Pilyug m & Thuan (2005)| l the authors re- 
vised the P-index calibration proposing two 
different functional forms for the upper branch 
and lower branch in the plane R23 vs log 0/H 
respectively. For the upper branch: 

12 + logo/// = 

(R23 + 726 A + 842.2 ■ P + 337.5 • P^) 
(85.96 + 82.76 ■ P + 43.98 • + 1.793 ■ R23) 

These relations can be used to determine the 
Oxygen abundances of the HII regions in SF 
galaxies and remove the degeneracies, since 
the P-index increases from left to right in the 

figure nn 

4.3.1 . Our sample 

Applying this procedure, we have found that 
our sample of SFs galaxies has P values in 
the range [0.1-0.9], consistent with the up- 
per branch sequence and that the sources are 
in the transition zone [12] At the same time, 
the estimated Oxygen abundance of our sam- 
ple is in the range [7.9-8.4]. Assuming the 
oxygen solar abundance of - 4.5710""^ 
( Asplund et al. (2005)] ), the gas in this sam- 
ple of galaxies has a sub-solar abundance. If 
the metal abundances scale with the oxygen 
abundance, how it is usually assumed, the gas 



Rafanelli: The nature of gas and stars in the circumnuclear regions of AGN 



297 




i^-«T lirnrnr. 



Fig. 11. Family of curves in the {R23, log 0/H) 
plane labeled by different values of the excita- 
tion parameter P, superimposed with data from 
different observations. The high-metallicity 
HIT regions with 0.0 < P < 0.3 are shown 
as opens squares, those with 0.3 < P < 0.6 
as plus sign and those with 0.6 < P < 0.9 as 
open circles. The low-metallicity HII regions 
with 0.5 < P < 0.7 are shown as filled trian- 
gles, those with 0.7 < P < 0.9 as crosses and 
those with 0.9 < P < 1.0 as filled circles. 

metallicity is between 0.2 and 0.5 Zq in fig- 
ure [T2] It needs to be stressed that these de- 
terminations of the metal abundances take into 
account only the gas component and that, in 
order to achieve the real metallicity, dust com- 
position has to be carefully modeled and con- 
sidered. Assuming that the ISM has solar abun- 
dance, about 40% of the oxygen resides in sil- 
icate grains, implying that = 1.9410"'* 
oxygen atoms per hydrogen atom are in dust 
form ( [Groves et al. (2006)l l. Even if there is an 
increase in the oxygen abundance, the metal- 
licity is in any case sub-solar or solar (the esti- 
mated value represents an upper limit because 
the dust content, usually, increases with the 
metallicity). 

5. Metallicity from stellar population 
synthesis 

The global relationship between the stel- 
lar populations and structural proper- 
ties of hot galaxies has been studied in 




Fig. 12. The star-forming galaxies of the sam- 
ple used in this work in the (7?23, log (9///). 
These galaxies have excitation parameter P in 
the interval [0.1,0.9], values that are consis- 
tent with their location in the upper branch se- 
quence of the plot. 




t 



Fig. 13. Same plot as before, with the assump- 
tion that total metallicity scales linearly with 
the Oxygen abundance. In this case, the SFs 
galaxies of our sample have metallicity be- 
tween 0.2 ■ and 0.5 ■ Zq 



( [Bender et al. (1993)l l using two different mea- 
sures of the global stellar population, namely 
the Mg2 index and the {B — V)q colour. The 
Mg2 index, introduced in ( [Faber et al. (1977)| ), 
is defined as the the diff^erence in mag- 
nitude between the instrumental flux in a 
window F2 = [5 156.0, 5 197.25] A centered 
on the Magnesium spectral feature and the 
pseudo-continuum interpolated from two 
windows. Fx = [4897.0, 4958.25] A and 
F3 = [5303.0,5366.75]Aatflieblueandred 
sides respectively: 



Mg2 = 



Fx +0.61(^3-/^1) 



(24) 



Using population synthesis models of var- 
ious authors, the relation between the Mg2 in- 
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dex and the central velocity dispersion cry can 
be translated in an approximate relation be- 
tween the same spectroscopic index and the 
age and metallicity of the stellar population: 

log Mg2 = 0.4 1 ■ log — + 0.4 1 ■ log r - 1 .00(25) 

Zq 

where Z is the metallicity, t is the age mea- 
sured in Gyr and the uncertainty on the coef- 
ficients is about +20%. This relation is rela- 
tively reliable for -0.5 < log (Z/Z©) < 0.3 
and 5 < f < 15 (always expressed in Gyr). 
The age and metallicities of the stellar pop- 
ulations of the normal galaxies considered in 
our work have been determined by exploiting 
this relation. In turn, one of the two param- 
eters has been derived from the fitted model 
of the continuum component of the spectrum 
obtained by STARLIGHT (whose output con- 
sist in a set of weights expressing the relative 
contribution of each single stellar population 
of given age and Z to the observed spectrum); 
this estimate of either Z or f has been used to 
calculate the other parameter as an unknown 
term in the |25] relation, using as approxima- 
tion of the measured values of the Mg2 values 
the equivalent width of the Mgl lines, retrieved 
from the SDSS database. Both temperature and 
metallicity of the stellar population could have 
been, in principle, derived by the only spectral 
models obtained by STARLIGHT, but this ap- 
proach may prove to be essentially wrong due 
to the intrinsic limited coverage of the (f, Z) 
plane of the library of template stellar popula- 
tion used by STARLIGHT to fit the continuum 
component of the observed spectra, since ev- 
ery possible final estimates of both parameters 
is a weighted average of age and metallicities 
of the stellar population contained in the tem- 
plate Ubrary. 

The age t and metallicity Z distributions 
for all galaxies belonging to the three differ- 
ent samples considered in this work and de- 
rived by STARLIGHT only are shown in fig- 
ures [14] and [15] as "whisker and box" plots 
with, on a side, the histograms obtained from 
the same distributions. As expected, the normal 
galaxies have significantly older stellar popu- 
lations than the other two samples of galax- 
ies, among which SFs are still younger than 



Age3 (liglit-Traetian weighted) 
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Fig. 14. Distribution of mean weighted age t 
for the three samples of galaxies considered in 
this work, obtained using STARLIGHT. The 
"box and whisker" plots show the position of 
the first four quartiles of the age distribution 
and the outliers. 




Fig. 15. "Box and whisker" plots of the distri- 
butions of the mean weighted metallicity Z for 
the three samples of galaxies considered in this 
work. 



Seyfert 2 galaxies. The metallicity distribu- 
tions show that normal galaxies and Seyfert 
galaxies, while almost consistent from a sta- 
tistical point of view, possess more evolved 
stellar populations than star-forming galaxies. 
On the other hand, the same distributions ob- 
tained using the measured values of the Mg2 
index and the relation [25] are shown in fig- 
ures ?? and [T7] respectively. The comparison 
between these plots indicates that, though the 
relative relations between distributions of dif- 
ferent samples of galaxies remain unchanged, 
stellar populations appear to be systematically 
younger and less metal-rich than in the case of 
the estimation by STARLIGHT only. 
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Fig. 16. "Box and whisker" plots of the distri- 
butions of the mean weighted age t for the three 
samples of galaxies considered in this work. 
These estimates have been obtained using the 
relation |25] the mean metallicity provided by 
STARLIGHT and the observed Mg2 indexes 
of the galaxies. 

Z Trorn Mgl (light -fraction weighted) 




Fig. 17. "Box and whisker" plots of the distri- 
butions of the mean weighted metallicities Z 
for the three samples of galaxies considered in 
this work. These estimates have been obtained 
using the relation |25j the mean weighted age 
provided by STARLIGHT and the observed 
Mg2 indexes of the galaxies. 

6. Conclusions 

In this work the metal content of the stellar and 
gaseous components of the central regions of 
galaxies has been characterized with multiple 
methods. Three different samples of galaxies, 
namely normal galaxies, star-forming galax- 
ies and galaxies hosting an AGN (Seyfert 2), 
retrieved from the SDSS DR7 spectroscopic 
database and classified according to an origi- 
nal method based on spectroscopic diagnostics 
(see |Vaona (2 009) ), have been compared. For 
what the nuclear stellar populations are con- 
cerned, the presence of a strong blue contin- 



uum and the absence of any 4000A break in 
the spectra of the SF galaxies are an indication 
of the presence, in their circumnuclear regions, 
of an excess of young hot blue stars compared 
with NG and S2G. On the contrary, the spectra 
of the circumnuclear regions of S2G and NG 
indicate quite similar conditions: a dominat- 
ing stellar component of type F and a 4000A 
break in NG deeper than in S2G of a factor 
~L5. This indicates that the presence of an ex- 
cess of A stars in the central regions of S2G 
compared to NG is quite likely, suggesting that 
there have been in the recent past of S2 AGN 
slightly more star formation events than in NG. 
The analysis of the gas of the central regions of 
the Seyfert 2 galaxies has been performed us- 
ing a very robust modeling of the photoioniza- 
tion processes occurring in gas clouds around 
the AGN and exploring a very large region 
of the model parameters space in order to re- 
produce realistically the observed spectra from 
simulation. Both single and two clouds models 
indicate for S2 a significantly larger metallicity 
than solar (L5 ^ 2.5 Z/Z© with SC models and 
1.0 2.0 Z/Zq for the 2C models). These re- 
sults support the idea that the AGN fuel is com- 
posed by evolved material and that the gas has 
a local origin in most of the observed Seyfert 
2 galaxies. The gas metal abundances for SF 
galaxies have been estimated using the method 
of the P-index (Pilyugin & Thuan (2005)) and 
we have obtained results which clearly in- 
dicate sub-solar metal abundances (0.2 - 0.5 
Zq), confirming that young stars are fueled 
by relatively primordial and unprocessed ma- 
terial. The overall relations between the av- 
erage age and metallicity of central stellar 
populations of the three different samples of 
galaxies here considered have been deter- 
mined using the Mg2 absorption feature as 
proposed by( |Faber et al. (1977)) and exploit- 
ing the spectral continuum components of ob- 
served spectra of the galaxies extracted using 
STARLIGHT ( jRafanelh et al. (2009)| l. This 
analysis qualitatively confirms that Seyfert 2 
galaxies have young low-metallicity circumnu- 
clear stellar populations feeding in a metal-rich 
environment where gas clouds are composed 
of rather high metallicity material; on the other 
hand, SF galaxies have very young and primor- 
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dial central stellar populations. In conclusion, 
the comparative study of metallic content of 
stars and gas in this work indicate that Seyfert 
2 nuclei occur preferentially in a stellar envi- 
rormient composed of young stars where recent 
circumnuclear star formation processes were 
active (as witnessed by the high metallicity of 
gas clouds), even if there is no compelling ev- 
idence that would connect present or past star 
formation activity to the nature of the central 
engine. The two phenomena, SF and AGN ac- 
tivity, may therefore be not causally connected 
to each other or, if so, through indirect or/and 
not clear physical mechanisms. 
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